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Abstract 

In this paper we calculate the quark number susceptibility (QNS) of QCD at finite temperature 
under the rainbow-ladder and Ball-Chiu type truncation schemes of the Dyson-Schwinger approach. 
It is found that the difference between the result of the rainbow-ladder truncation and that of Ball- 
Chiu type truncation is small, which shows that the dressing effect of the quark-gluon vertex on 
the QNS at finite temperature is small. It is also found that at low temperature the quark number 
susceptibility is nearly zero and it increases sharply when the temperature approaches the chiral 
phase transition point. A comparison between the result in the present paper with those in the 
literature is made. 
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I. INTRODUCTION 



As is well known, dynamical chiral symmetry breaking (DCSB) and confinement are 
two fundamental features of quantum chromodynamics (QCD). Nowadays it is generally 
believed that at high enough temperature and/or density strongly interacting matter will 
undergo a chiral symmetry restoring and deconfining phase transition (for a recent review, 
see for example Ref. The study of such phase transitions has become a field of great 
theoretical and experimental interest over past years because it will provide us a fundamental 
understanding of the basic theory including the origin of observed mass and the nature of 
the early universe. A central goal of relativistic heavy ion collision experiments is to explore 
the phase structure and investigate the chiral phase transition of QCD matter. 

In the study of chiral phase transition the enhanced quark number fluctuations are 
thought to be an essential characteristic In the confined/chirally broken phase 

quark (baryon) numbers are associated with hadrons in integer units whereas in the de- 
confined/chirally restored phase they are associated with quarks in fractional units which 
could lead to different quark number fluctuations in the two phases. Theoretically such fluc- 
tuation is constructed from measurement of quark number susceptibility (QNS) as a function 
of temperature, which is the the response of qusirk number density to infinitesimal change of 
chemical potential at zero chemical potential |2|]. In particular, it was recently argued that 
the QNS may be used to identify the chiral critical point in the QCD phase diagram js-lOj. 



Hence the QNS of QCD has been extensively studied over the past twenty years using dif- 



ferent approaches including lattice QCD simulations 



E 



11 



16 



, Nambu-Jona-Lasinio model 



5|, |l7|, hard-thermal/dense-loop resummation techniques |18|-|22| and the rainbow-ladder 



approximation of the Dyson- Schwinger equations (DSEs) approach |23|-|25| . It should be 
noted that many model calculations of the QNS are done at the level of mean field approxi- 
mation (it corresponds to the rainbow-ladder approximation in the framework of the DSEs 
approach) 23l-l25l] . Since the QNS is a very important quantity in characterizing the chiral 
phase transition, theoretically it is very interesting to calculate this quantity beyond mean 
field approximation, i.e, beyond the rainbow-ladder approximation of the DSEs approach. 
The purpose of the present work is to calculate the QNS at finite temperature beyond the 
rainbow-ladder approximation of the DSEs approach. 



Over the past few years considerable progress has been made in the framework of the 
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rainbow-ladder approximation of the DSEs approach 26l-l3l[| (for a recent review, see Ref. 



32|). Due to the great success of the rainbow-ladde r ap proximation (RL) of the DSEs ap 



proach in hadron physics, the authors in Refs. [23|, [2^ adopt this approximation to solve 
DSEs for the dressed quark propagator and the vector vertex and from these obtain the nu- 
merical value of the QNS. However, it is well known that the rainbow-ladder approximation 
uses a bare quark-gluon vertex, which violates Slavnov- Taylor identity of QCD. In order to 
overcome this deficiency, physicists are trying their best to go beyond the rainbow-ladder 
approximation. Much work has been done in this direction, including the Ball-Chiu (BC 



33 



3^, the CP vertex [35 1 



vertex derived from the vector Ward-Takahashi identity (WTI) 
which takes into account some transverse effects and the vertex derived from the transverse 
WTI 36|-|38|. As was shown in Ref. [39|], if the ghost amplitudes and the gluon dressing 



function factor are ignored in Slavnov- Taylor identity one would find the resultant relation 
has the form of a color matrix times the WTI structure which could be satisfied by the 
BC vertex ansatz multiplied by the color matrix. Therefore, in this paper we adopt such 
a quark-gluon vertex ansatz to explore the effect of vertex dressing on the quark number 
susceptibility. 

When the BC vertex is used to calculate the dressed quark propagator from the DSEs, 
one should construct a consistent kernel approximation corresponding to this vertex. How to 
construct systematic and convergent expansions for the kernels of DSEs is a long-standing 
unsolved problem. Recently, an important progress in this aspect has been achieved in 
Ref. j4o|. The authors in Ref. 40 1 have proposed a Bethe-Salpeter kernel which is valid 
for a general quark-gluon vertex and this provides a theoretical foundation for calculating 
the QNS beyond the rainbow-ladder approximation. In the present paper we will use this 
method to calculate the QNS at finite temperature. 



II. GENERAL FORMULA FOR QNS 

To make this paper self-contained, let us first recall the definition of the QNS which is 
analogous to the familiar magnetic susceptibility. If we confine ourselves to the two-flavor 
case with exact isospin symmetry and set /i« = /id = [fiu and fid are the chemical potential 
of the up and down quarks), the quark number density and the corresponding susceptibility 
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are given by 



and 



X{T) 



T din Z{fx,T) 
dpifi, T) 



dp 



(1) 



(2) 



^l=o 



respectively, where Z{p,T) is the partition function of QCD at finite temperature T and 
quark chemical potential p. 

From Eq. ([T]) and by using functional integral techniques, one can derive a well-known 
result (see, for example, Refs. |4ll. |42|) 

p{p,T) = -N^NfT J 7A-,iHG[p]{Pkh,} , (3) 

where G[p]{pk) is the dressed quark propagator at finite T and p, pk = {p,uJk) and = 
{2k + 1)ttT are the fermion Matsubara frequencies with k being integers, Nc and Nf denote 
the number of colors and of fiavors, respectively, and the trace operation is over Dirac indice. 
From Eq. ([3]) it can be seen that the quark number density p{p, T) is totally determined by 
the dressed quark propagator at finite chemical potential and temperature. If one substitutes 
the free quark propagator at finite T and p into Eq. ([3]), one will obtain exactly the Fermi 
statistics result for the quark number density of a free quark gas, as was shown in Ref. 42|. 
Substituting Eq. ([3]) into Eq. ([2]), one finds the following expression 

d^p jdG[p]{pk) 

By means of the following identity 

dG-'[p]{pk) 

^l=o 



X(T) 



dp 



-74 



(4) 



/i=0 



dp 



-r4(Pfc,o), 



(5) 



where T^ip^, 0) is the fourth component of the dressed vector vertex r^(j9fc, 0) (the chemical 
potential p can be regarded as the fourth component of a constant external vector field 
which is coupled to the vector current of quarks, see Ref. |23|), one can find the following 
general formula for the QNS at finite T 



X{T) 



r d^n 

NcNfT J j^tT[G{p,Uk)T,{p,,0)G{p,u,h,], 



(6) 



which is the same as the expression given in Ref. 23|. It should be noted that in contrast 



to the chiral susceptibility 



431], there is no ultraviolet divergence in the quark number 



susceptibility, and if one substitutes the free quark propagator (in the chiral hmit) at finite 
T and the bare vertex into Eq. ([6]), one will obtain exactly the QNS of a free massless quark 



;as Xfree = NfT^ |23|. 



Since Eq. (jH]) is a model-independent expression for the QNS at finite T, one can obtain 
the exact result of the QNS at finite T once the exact dressed quark propagator and vector 
vertex at finite T is known. However, at present it is very difficult to calculate the dressed 
quark propagator and vector vertex from first principles of QCD and hence one has to resort 
to various nonperturbative QCD models. In the following we will use DSEs of QCD to 
calculate the QNS from Eq. 

III. DSES APPROACH 

At finite temperature T the Dyson-Schwinger equation for the quark propagator can be 
written as follows (the renormalization constants are set to one and this will be explained 
later) 

where Gq^ is the inverse of the free quark propagator, g is the strong coupling constant, 
Df^i, is the dressed gluon propagator, is the dressed quark-gluon vertex, g„ = (g, Un) and 
Un = (2n + l)7rT with n being integers. At finite temperature T the 0(4) symmetry of QCD 
is broken to 0(3) (in the present paper we shall always work in Euclidean space). From 
general Lorentz structure analysis one finds that the inverse of the dressed quark propagator 



G-\pk) = G,\p,) + \TY.j -^,g'D^,{p, - qn)l,G{<ln)n{<ln.Vk), (7) 



at finite T can be decomposed as follows 
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G-\pk) = ii-pA{f.^l)+^^kl,G{f,ul) + B{f,ul), (8) 

where A, B and G are scalar functions of fP' and w^. 

The dressed vector vertex V ^ satisfy the following inhomogeneous Bethe-Salpeter equation 



TM^Pnd = + J ^[0(g+)r^(g„,PnJO(g_)],,ir[J(g„,pfc;PnJ, (9) 

— * 

where pk is the relative and Pq, the total momentum of quark-antiquark pair, Pq^ = (P, Qi) 
and Qi = 21'kT with / being integers, g± = g„ ± PnJ'i and r, s, t, u represent color and Dirac 
indices. Here K is the fully-amputated quark-antiquark scattering kernel. According to 



Eq. (|6]), when calculating the QNS we need to know the fourth component of the dressed 
vector vertex at zero total momentum T4{pk,0) which can be decomposed as follows (see 
the appendix): 

r4(Pfc, 0) = a;fc(Fi7 ■ p + F^ujklA - iF^), (10) 

where Fi, F2 and -F3 are scalar functions of fP' and w^. For the free field case, Fi = F3 = and 
F2 = l/ul Substituting Eqs. (jH]) and into Eq. (ED, one finds the following expression 
for the QNS at finite temperature 

2NJ^ foo , , 2BCF, + 2ACF^f + F2iC'ul-A'f-B^) 

and if one uses the free field value A = C = 1, B = Q and Fi = F3 = 0, F2 = l/w^ in Eq. 
ffTTj) one will re-obtain the free quark gas result. 

When one tries to solve Eqs. ([7]) and ([9]), one needs to input the model gluon propagator 
in advance. In the present work we employ the following model gluon propagator at finite 
temperature {Q^i = {Q,^i)): 

9'D,,{Qn,) = Pl,{Qn,)DT + PiAQn,)DL, (12) 

where -Pj^^ are the transverse and longitudinal projection operators, respectively: 

_ . _ Qi^Qu , . 



fj, and/or z/ = 4 



Pl= { ,-2 ' (14) 

5ij-QiQjlQ^ /i, z/ = 2, j = 1,2,3 



and 



Dr = l^;^Q^^,e-«.'"^ (16) 
fl. = ^^(«,+"J)e-"'°. (17) 
This model gluon propagator is a simplified form of the effective interaction proposed in Ref. 



451], which is the finite temperature extension of Maris- Tandy model 26|, \2j\. It should be 
noted that this model delivers an ultraviolet finite gap equation for the quark propagator 
and hence the regularization mass-scale could be removed to infinity and the renormalization 



constants set to one. In this model = IGvr^T^/S is a temperature-dependent mass-scale 
and uj and D are active parameters of the model which are not independent: a change in D 
will be compensated by an alternation of u [29]. For uj G [0.3, 0.5] GeV in the rainbow-ladder 
truncation scheme, fitted in-vacuum low-energy observables are approximately constant if 
ujD = (O.SGeV)'^ and in the present paper we use u = 0.5 GeV. 

Now let us discuss the truncation schemes of DSEs. Under the rainbow approximation 
one uses bare vertex 'ji, to replace and Eq. ([7]) becomes (in the chiral limit) 



4 r d'^q 



(2v 



(18) 



Under the corresponding ladder approximation for the fully-amputated quark-antiquark 
scattering kernel K one uses 



A" 



(19) 



tT 



and Eq. ([9]) becomes 

F^(pfc,PnJ = 7/.- 3^1] y (27r)3^ Dp„{pk- qn)lpG{q^)V^{qn,Pa^)G[(lJ)^„. (20) 

The rainbow-ladder approximation is the lowest order truncation scheme for the DSEs 46| . 
and due to the reasons mentioned in the introduction physicists are trying to go beyond it 
for years. Here the key points are the dressed quark-gluon vertex and the corresponding 
quark-antiquark scattering kernel. Just as is shown in the introduction, in this work, for the 
dressed quark-gluon vertex we will employ the following finite temperature extension of 



BC vertex 



33 



34| 



(r 



BC -rBC'\ 

5 4 ) 



where 



f^^(g„,Pfc) = EA7 + {q + p) 



7 • (g + P) ^ , 74(Wn + Wfc) 

Aa H Ac - «Ab 



with ( = A, B, C) 



Sc 74 + (w„ + Uk] 



1-{q + P) ^ , 74K + CJfc) 
H — 
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(21) 

(22) 
(23) 

(24) 
(25) 



Therefore Eq. ([7|) becomes 

G^^iPk) = ii ■p + iujkl4 + :^TJ2 J j^^9'^D^y{pk - qn)l^G{qn)T^^ {qn,Pk)- (26) 

Now one should find a kernel for Eq. ([9]) which is consistent with BC vertex. This is a 
difficult task and recently great progress has been done on this aspect. The authors in Ref. 

find a way to constrain the kernel for a general vertex. For BC vertex, following their 
method, Eq. ^ can be written as 

4T f d?q 

+ y J^;^9^Dp^iPk - qnhpG{qn)Af,a{Pk, qn, 0), (27) 

where A^,n-(pfr, qr?,; 0) is a four-point Schwinger function which is completely defined by the 



quark self-energy 48|, |49|. It satisfies a similar identity as those in Ref. 40 1 



iiPk - qn)a■A^,a{Pk, qn] 0) = T^{pk, 0) - r^(g„, 0), (28) 

and from this identity one can write down the expressions of A^o-(pfc, qn] 0) (for details, see the 
appendix). Now, by means of the BC vertex ansatz and the corresponding quark-antiquark 
scattering kernel in Eq. fl27|) (in this paper we shall call this truncation scheme a BC-type 
truncation) , we can numerically calculate the dressed quark propagator and the vector vertex 
at finite T beyond the rainbow-ladder approximation, which is needed to calculate the QNS 
of QCD at finite temperature. 



IV. NUMERICAL RESULTS 



By using numerical iteration method one can solve Eqs. ( fT8l) . ( I20l) . ( |26|) and ( 1271) . and 
the results are shown in Figs 1-4. Here it should be noted that when working with the 
rainbow- ladder truncation, we employ model parameters in which D = 1.0 GeV^. Now 
when calculating the dressed quark propagator using the BC-type truncation, we employ 
the model gluon propagator of the same form. However, because the amount of chiral 
symmetry breaking (as measured by the chiral condensate) and related quantities such as 
the pion decay constant are very different between the rainbow-ladder truncation and the 
BC-type truncation [40], one should employ refitted model parameters in Eqs. (fT6l) and ( 1T7|1 
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50 



5l| when performing the calculation in the BC-type truncation. Under the BC vertex, 



;he va. 



50 



ue of the parameter D fitted from the chiral condensate is D = 0.5 GeV (see, Refs. 



51|). and in this paper we choose this value. 
The dressing functions A{0, tc^T"^), C(0, tt^T^) and -8(0, tc'^T'^) as functions of temperature 
are shown in Fig. 1 and Fig. 2, respectively. If one takes -B(0,7r^T^) as the chiral order 
parameter, from Fig. 2 it can be seen that the transition temperature Tc ~ 102 MeV for 
both RL and BC-type truncations, because at this temperature -6(0, vr^T^) for both two 
cases decrease to zero. This temperature is smaller than the result of lattice QCD in which 
Tc ~ 150 MeV (see, for example, Ref. 52|) and this discrepancy may be ascribed to the 
fact that we have ignored the perturbative tail in our model gluon propagator (IT^ and (ITTI) . 
From Fig. 1 it can be seen that the derivative of A{0, tc^T'^) and C(0, tt^T^) with respect to 
T undergo a sudden change at Tc. The temperature dependence of -Fi(0, tt^T"^), -^2(0, vr^T^) 
and -^3(0, vr^T^) is shown in Figs. 3-4. One can see when T > T^, the main dressing effect of 
the vector vertex T4{pk,0) comes from F2. This is reasonable because from Eq. ( fTOl) 
is the coefficient of the bare vertex 74, which should be a dominant term when T is high 
enough. 

The calculated QNS under both truncation schemes are shown in Fig. 5 and Fig. 6. In 
Fig. 5 xbc and xrl are the QNS obtained using BC-type truncation and RL truncation, 
respectively. {qq)BC and {qq)BCO are the chiral condensate obtained using BC-type trun- 
cation at finite T and zero T, respectively. {qq)RL and {qq)RLo are the chiral condensate 
obtained using RL truncation at finite T and zero T, respectively. It can be seen from Fig. 
5 that when T < Tc, the QNS increases with increasing T and the rate of increasing also 
increases with T. At T = Tc, x/Xfree ~ 0.95 and the derivative of QNS with respect to T 
undergoes a sudden change (from positive value to negative value). Such a behavior could 
be regarded as the signal of chiral phase transition. When Tc < T < 200 MeV, the QNS 
decreases slowly with increasing T. From Fig. 6 it can be seen that when T is larger than 
about 200 MeV, the QNS again increases very slowly towards the free quark gas limit with 
increasing T, which is the result of asymptotic freedom of QCD. In all ranges of the tem- 
perature investigated, the QNS obtained using BC-type truncation is slightly larger than 
the one obtained using RL truncation, but the difference is smaller than 10% when T > 95 
MeV, which indicates that RL approximation is good for calculating the QNS at finite T. 
The chiral condensates obtained using both truncation schemes are also shown in Fig. 5, 



9 



2.2 








2.1 


. A(0,7iT')ofBC 






2.0 


C(0,tiTVBC 






1,9 








1,3 


C(0,Tt'T'')of RL 






1,7 








1.6 
1.5 








1.4 








1.3 








1.2 








1.1 









0.12 0. 
T (GeV) 



FIG. 1: A{0, vr^T^) and C(0, vr^T^) of the 
quark propagator. 




FIG. 2: B{0,tt'^T^) of the quark propaga- 
tor. 
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FIG. 3: Fi(0,7r2r2) and F2 (0,^2^-2) 
the vector vertex. 
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FIG. 4: F3{0,iT^T^) of the vector vertex. 



in which they are normahzed, respectively, by the values of the zero temperature results 
obtained using the corresponding truncation schemes. From Fig. 5 it can be seen that 
when T < Tc the chiral condensate decreases monotonously with increasing temperature, 
and when T > Tc it equals zero. One can also see that the chiral condensates obtained using 
the two truncation schemes are almost the same, which shows that the dressing effect of the 
quark-gluon vertex on the chiral condensate is very small. 

It is instructive to compare our result of QNS with the corresponding results given in 



Refs. 2^ and 25|. In Ref. 25|] a parametrized quark propagator which is based on an 
effective interaction similar with the one used in the present paper is employed and the 
resultant QNS exhibits similar behaviors with the result of the present paper. However, the 
QNS given in Ref. [25| becomes negative when 50 MeV < T < 90 MeV, and according to 
the result of the present paper, such a behavior might be an artifact of the model quark 
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FIG. 5: Temperature dependence of the 
QNS and the chiral condensate. 



FIG. 6: The QNS in the range T > 200 
MeV 



propagator used in Ref. 



25|. On the other hand, in Ref. 2J] a rank-1 separable model for 



the effective gluon propagator is employed and the resultant QNS has a peak at Tc. When T 
is across the derivative of the QNS with respect to T changes sign, which is in agreement 
with the result in the present paper. However, when T is around T^, the QNS calculated in 
Ref. is larger than the free quark gas result. This point is hard to understand physically, 
because the free quark gas result should be an upper limit and a result beyond this limit 
means the space-like damping mode is extremely dominant in the energy spectrum of quark 
matter. According to the result of the present paper, such a behavior might be an artifact 
of the model gluon propagator choosed in Ref. 24] • addition, from Fig. 5 it can also 
be seen that for both RL and BC trancation schemes, the bahavior of the QNS across the 
phase transition point is nonmonotonic. Such a nonmonotonic behavior of the QNS can be 
regarded as a prediction of our model. Here it should be pointed out that the nonmonotonic 
behavior of QNS around Tc obtained in our model differs from the existing lattice results 



of Refs. 



15 



161] where a monotonic behavior of the QNS across the phase transition point 



is found. It should be noted that in the present work we are working in the chiral limit. 



whereas in the lattice simulations of Refs. 



15 



16| a nonzero current quark mass is assumed. 



It is likely that the nonmonotonic behavior of QNS around in the case of chiral limit 
obtained in the present work is not within reach of the existing lattice studies. 
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V. SUMMARY 



In summary, we calculate the QNS of QCD at finite temperature using both RL and 
BC-type truncation schemes of DSEs approach. Our result shows that when T < T^, the 
QNS increases with increasing T and the rate of increasing also increases with T. At 
the derivative of QNS with respect to T undergoes a sudden change (from positive value to 
negative value). Such a behavior can be regarded as the signal of chiral phase transition. 
When Tc < T ^ 200 MeV, the QNS decreases slowly with increasing T. When T is larger 
than about 200 MeV, the QNS again increases very slowly towards the free quark gas limit 
with increasing T, which is the result of asymptotic freedom of QCD. When T ^ 95 MeV, 
the difference between the QNS obtained using RL truncation and the one obtained using 
BC-type truncation is small (< 10%), which shows that the dressing effect of the quark- 
gluon vertex on the QNS is small. In addition, from our study it is found that for both 
RL and BC truncation schemes, the bahavior of the QNS across the phase transition point 
is nonmonotonic. Such a nonmonotonic behavior of QNS around Tc obtained in our model 



differs from the existing lattice results of Refs. [15|, Il6| where a monotonic behavior of 
the QNS across the phase transition point is found. This difference deserves further study. 
Finally, we want to stress that the method used in the present paper can be generalized to 
the case of both finite temperature and finite chemical potential and be used to investigate 
the properties of the critical end point and the phase diagram of QCD. 
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Appendix A: Expressions of the BC-type kernel 



From general Lorentz analysis one can find the vector vertex r^(pfc, 0) has the following 
structure (z = 1, 2, 3) 



= Wfc(Fi7 ■ p + F2UJklA - iFs) 



(Al) 
(A2) 
(A3) 



where F2 = F2 + Fq/uI, Eq, . . . ,E3 and Fq, Fi, F2, F2, F3 are scalar functions of and ul- 
For the free field case, Eq = Fq = 1, Ei = E2 = E^ = Fi = F2 = F^, = and hence 
F2 = 

From Eq. (1251) the four-point Schwinger function A^^-iPk, In] 0) can be written as follows 
(^,J = 1,2,3) 



iAij{pk,qn;0) 



- i^EsiPi + Qi) 



2AiJo7i + {pa ■p + qii-q) + ^E^iPiUJk + qi(^nh4 

Pi + Qi , 7 • (P + ^) X 



El 



iAi4(pfc,g„;0) 



iAii{pk,qn;0) 



iA44(pfc,gn;0) 



+ ^E,- 

^k + 



2AEoji + Ae^ {pa ■p + qii-q) + AE.,{PiUJk + qi(^n)iA 



iAEsiPi + qi)] + T,E2 ^' 't ^' 74 



Pi + qi 



2 

- iAp.^ (uk + uj, 

^k + 



(c^fc7 ■ P + ^nl ■q) + Ap^ {ujI + W^)74 



-li 



(A4) 



(A5) 



(A6) 



A^i (Wfc7 ■ P + ^^n7 ■q)+ (Wfc + W^)74 

7 • + 



iAF.,{uJk + C^n) 



Si^^ (wfc + w„)74 - iSi^'g (A7) 
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